Dependence of the electronic structure of self-assembled InGaAs/GaAs quantum dots 

on height and composition 
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While electronic and spectroscopic properties of self-assembled Ini_ x Ga x As/GaAs dots depend on 
their shape, height and alloy compositions, these characteristics are often not known accurately from 
experiment. This creates a difficulty in comparing measured electronic and spectroscopic properties 
with calculated ones. Since simplified theoretical models (effective mass, k-p, parabolic models) 
do not fully convey the effects of shape, size and composition on the electronic and spectroscopic 
properties, we offer to bridge the gap by providing accurately calculated results as a function of the 
dot height and composition. Prominent features of our results are the following, (i) Regardless of 
height and composition, the confined electron energy levels form shells of nearly degenerate states 
with a predominant orbital "s", "p", • • • character. In contrast, the confined hole energy levels form 
shells only in flat dots and near the highest hole level (HOMO), (ii) In alloy dots, the electrons' 
"s-p" splitting depends weakly on height, while the "p-p" splitting depends non-monotonically — due 
to alloy fluctuations. In pure, non-alloyed InAs/GaAs dots, both these splittings depend weakly on 
height. Further, the "s-p" splitting is larger while the "p-p" has nearly the same magnitude. For 
holes levels in alloy dots, the "s-p" splitting decreases with increasing height (the splitting in tall dots 
being about 4 times smaller than in flat dots), whereas the "p-p" splitting remains nearly unchaged. 
Shallow, pure non-alloyed dots have a "s-p" splitting of nearly the same magnitude, whereas the 
"p-p" splitting is about three times larger, (iii) As height increases, the "s" and "p" character of 
the wavefunction of the HOMO becomes mixed, and so does the heavy- and light-hole character, 
(iv) In alloy dots, regardless of height, the wavefunction of low- lying hole states are localized inside 
the dot. Remarkably, in non-alloyed InAs/GaAs dots these states become localized at the interface 
as height increases. The localized states are nearly degenerate and polarized along [110] and [110]. 
This localization is driven by the peculiarities of the biaxial strain present in the nanostructure. 



I. INTRODUCTION 



The electronic structure and spectroscopic prop- 
erties of quantum dots, including excitonsjiiSiSi**** 
charged excitons^i^ multiexcitons^i^ and excitonic fine- 
structure^ all depend on the size and shape of the dots. 
This dependence reflects both quantum confinement ef- 
fects, as well as shape-induced band-folding and inter- 
band coupling. 7 Simulations of electronic structure and 
spectroscopic properties of quantum dot a 8 i 9 i 10 i 1:L i 12 i 13 i 14 
must naturally assume the size and shape of the dot. On 
the other hand, measurements of spectroscopic properties 
of a ^0^2*^^^ are rarely accompanied by accurate mea- 
surements of the size and shape of the GaAs-covered dot, 
except in rare cases where detailed cross-sectional scan- 
ning tunneling micros copy experiments are performed, 
such as in Refs. Il5lll6lll7l This situation creates a sig- 
nificant difficulty, if not a crisis, in interpreting spectro- 
scopic data on quantum dots, and in critically testing 
various theoretical approaches. Thus, in reality one is 
often forced to address the inverse problem^ namely, fit 
the spectroscopic data to a theoretical model by using the 
size and shape as adjustable parameters. The difficulty 
with this approach is threefold: First, in this approach all 
theories, no matter how naive, ultimately work by virtue 
of forcing a fit to experiment, even if the assumptions 
entering the theory may seem unjustified in their own 
right (e.g. assuming single-band effective mass; neglect of 
strain; neglect of spin-orbit coupling.) Second, there are 



usually too many free-parameters, which involve not only 
non-trivial shapes, but also unknown composition profiles 
(e. g., lni_ x Ga x As dots). Third, since the relationship 
between shape and spectroscopic properties is model de- 
pendent, unrealistic shapes are often deduced. For exam- 
ple, in simple effective-mass ^iSiiSiiS k-p^S or parabolic 
models the "p" and "d" levels of electrons and holes 
are degenerate if one assumes spherical, lens-shaped, cu- 
bic, or cylindrical dots. In contrast, in more advanced 
atomistic models — like empirical pseudopotentialSLSiS^ 
or tight binding^ — those levels are split even for the 
above mentioned ideal shapes, resulting in clear spec- 
troscopic signatures. To fit measured spectroscopic sig- 
natures of actual dots by simple theoretical models, one 
needs to assume at the outset irregular shapes. For in- 
stance, Dekel et alH needed to assume a parallepipedal 
box to explain their multiexciton data on non-alloyed 
InAs/GaAs dots; and Ferreira assumed shape distortions 
to explain fine structured Such assumptions are not 
needed in atomistic approaches to modeling. 

In this work we have used a high-level atomistic ap- 
proach to predict the spectroscopic characteristics of 
Ini_ x Ga x As/GaAs self- assembled dots as a function of 
the most crucial geometric parameter, namely the heigh. 
We calculate strain profiles, "p-" and "d-" level splittings 
as well as electron and hole wavefunctions. While ulti- 
mately it will be necessary for experimentalists to report 
the size, shape, and composition profile to which their 
spectroscopic data correspond, the type of study reported 
here may be used to bridge, in the interim, spectroscopy 
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Sketch of a lens-shaped InGaAs/GaAs quantum dot 




FIG. 1: Sketch of the geometry of self-assembled 
Ini_ x Ga x As/GaAs quantum dots. A section of the GaAs ma- 
trix (grey) has been ommited for clarity. The dot (white) is 
lens shaped with base diameter b and height h. The dot lies 
on top of a 1-monolayer-thick Ini_ x Ga x As wetting layer. Two 
auxiliary features are shown: (i) A plane that bisects the dot 

(dashed lines) and (ii) the line OO' , which pierces the dot 
through the center along the [001] direction. 

with theory without clouding the issue by severe theoret- 
ical approximations. Similar s tudies were c arried out by 
the group of Bimberg in Refs. I2rj27l28l29l by Shumway 
et al. in Ref. EI_by Williamson et al. in Ref. by 
Kim et al. in Ref. |3(J and by Pryor in Ref. |^ 



II. METHODS AND THEIR ILLUSTRATION 

A. Choice of dot geometries 

Among several geometries, Stranski-Krastanow quan- 
tum dots growth in lens shape. 32 Hence, we focus on 
lens-shaped, self-assembled Ini_ x Ga x As/GaAs quantum 
dots (QDs). In addition, we include a 1-monolayer-thick 
Ini_ x Ga x As wetting layer (WL). Figure 1 shows a sketch 
of the geometry of the nanostructure. (QD+WL+GaAs 
matrix). We focus on dots with x = 0.4 and pure InAs. 
The QDs have circular base with diameter b — 252 A and 
height h in the interval 20 A-75 A. 



B. Strain relaxation done via atomistic, not 
continuum elasticity 

The position-dependent strain profiles present in the 
nanostructure are usually addressed within (harmonic) 
continuum elasticity (CE) theoryi2ii2i Notwithstanding, 
strain profiles can also be calculated within atomistic 
elasticity (AE) in the form of a valence force fields In 
AE the strain energy is expressed in terms of atomic po- 
sitions and interatomic force constants. While the force 
constants are fit to the elastic constants Cn, Cn and 



C44, much like in CE, expressing the strain in terms of 
atomic positions contains much more information than 
expressing strain in terms of CE coordinates. For ex- 
ample, a square-based pyramid has a C± v point group 
symmetry if described by CE. This means that the [110] 
and [110] facets are taken to be equivalent. In this case, 
all "p" energy levels are degenerate. However, in AE the 
point group symmetry of a square-based pyramid made 
of InAs is C^v In this case, the [110] and [110] facets are 
symmetry-inequivalent and "p" levels are split. 

In general, three interrelated sources are responsible 
for the splitting of the "p" states^ A quantum dot 
with a base that globally possesses inversion symmetry 
in the (001) plane — like a lens, a truncated pyramid or a 
truncated cone — has in reality the lower C^ v symmetry 
(for pure, non-alloyed InAs/GaAs dots), which originates 
from the underlying zinc-blende atomic lattice. Thus, 
the first source is the lack of inversion symmetry in the 
Ci v point group, which manifests itself around the dot- 
matrix interfaces where the [110] and [110] directions are 
inequivalent. The second source is the propagation of 
the strain field towards the center of the dot as a conse- 
quence of the atomic relaxation. The third source is the 
piezoelectric effect (which magnitude is presently under 
debate 3 ^) that arises from the strain field of Civ sym- 
metry. It should be noted that approximations used in 
previous calculations of the piezoelectric effec t 3 ^' 38 have 
been shown to be crude and further investigations have 
been called foriSi 

Our choice of AE is based on a generalization of Keat- 
ing model to 3 terms — bond-bending, bond stretching 
and their cross terms>i^ We fit the elastic constants Cn, 
C12 and C44 of zinc-blende InAs and GaAs, as well as the 
correct dependence of the Young's modulus with pres- 
sure for both materials. In the Ini_ x Ga x As alloy sys- 
tem, the bond-stretching and the bond-stretching/bond- 
bending cross-term parameters for the mixed cation Ga- 
As-In bond angle are taken as the algebraic average of 
the In-As-In and Ga-As-Ga values. While the ideal bond 
angle is 109° for the pure zinc-blende crystal, to satisfy 
Vegard's law for the alloy volume, the value of 110.5° was 
used for the mixed bond angle. AE is superior to CE 
in that it does not assume harmonicity — in fact, anhar- 
monic effects can be explicitly included 3 ^ into the valence 
force fields 

Here, we calculate the position-dependent strain tensor 
e(R) within the atomistic elasticity approach. To calcu- 
late the strain tensor cubic components e,-j = x, y, 
z), we proceed in two steps: (1) We relax the atomic po- 
sitions within the supercell in order to minimize the elas- 
tic energy, which is given by a generalized valence force 
fields (2) We relate the relaxed (equilibrium) atomic 
positions with the unrelaxed atomic positions via strain 
tensor^4 At each equilibrium position R/ of atom I we 
identify the tetrahedron formed by its 4 nearest neigh- 
bors. This tetrahedron is distorted in comparison to the 
unrelaxed tetrahedron. Thus, the three edges of these 
tetrahedra that are determined by the vectors connect- 
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ing the four neighbors can be related by the strain tensor e as follows. 
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where R a /, R&j, and R c ; are the 3 vectors that connect, 
respectively, neighbors 1 and 2, 2 and 3, and 3 and 4 
in the equilibrium, distorted tetrahedron that encloses 
atom I. R° ; , R^, and R° ; are the corresponding vectors 
(edges) in the unrelaxed tetrahedron. From Eq. QJ, 
the cubic strain tensor components are calculated by a 
matrix inversion. 

Figure [2 shows the biaxial strain 
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\J {£xx ~ SyyY + {£yy ~ S zz f + (e zz ~ S xx ) 2 (2) 

and the hydrostatic distortion (isotropic strain) I = 
Tr(e) as a function of dot height. The results presented 
correspond to the strain values on a planar section in 
the [100] direction that bisects the dot (see Fig. [[). 
The spikes in the dot and wetting layer region are a 
consequence of the random arrangement of In and Ga 
atoms. Thus, the spikes structure significantly changes 
from dot to dot. Two additional features are also promi- 
nent, (i) The biaxial strain at the dot-matrix interface, 
located near the top of the dot, increases with height 
and slowly decays off. This decay is faster as the dot 
becomes taller. The latter can be understood by con- 
sidering the dot as a spherical cap inscribed in a sphere 
of radius p = [h 2 + (b/2) 2 }/2h (the taller the dot, the 
smaller p ) and by assuming that the decay of B out- 
side the dot is qualitatively well described by (po/p) 3 , 
which is the radial dependence of the biaxial strain out- 
side a continuum, elastic sphere of radius po embedded in 
a lattice-mismatched medium^ Here, p is the radial dis- 
tance from the sphere surface. In our simulations, Tr(e) 
shows a similar trend as B. (ii) Due to the underlying 
zincblende atomic structure, the biaxial strain drops to 
zero along the [Oil] and [Oil] directions, as well as along 
their equivalent crystallographic directions. (See white 
arrows in Fig. [21 ) This feature is robust upon changes 
in height, being present for all the quantum dots we con- 
sidered. 



C. The single-particle electronic structure is 
calculated via pseudopotential plane-wave method 
instead of k ■ p 



approach combines a pseudopotential description of the 
single-particle Hamiltonian with the linear combination 
of bulk bands (LCBB) method to solve for the energies 
and wavefunctions^i In this method, the Hamiltonian 
reads 



V 



so 



J2 5> a (R-R?;g), (3) 



a— In,Ga,As I 



where (3 is an empirical parameter that accounts for non- 
locality effects; V so is a non-local empirical operator 
that describes the spin-orbit interaction^ v a (r);E) is a 
screened pseudopotential (for atom of type a) that de- 
pends on strain; and R" is the relaxed vector position of 
atom I of type a. The dependence of the atomic pseu- 
dopotential on strain transfers to the electronic Hamil- 
tonian the information on atomic displacements. The 
strain-dependent pseudopotential reads 



i; Q (R-Rr;0)[l+ 7 aTr(?)] ) 



(4) 



where 7 a is a fitting parameter. It should be noted 
that w Q (R — Rf;e) is fit to bulk properties of GaAs 
and InAs, including bulk band structures, experimental 
deformation potentials and effective masses, and LDA- 
determined band offsets. In order to improve the trans- 
ferability of pseudopotential w Q (R — R";e), a simple de- 
pendence on the chemical environment of atom a is in- 
troduced. For instance, for a = As in an environment of 
p Ga atoms and p — 4 In atoms we use 
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The pseudopotentials used in this work have been suc- 
cessfully tested for quantum wells 

The wavefunction of state i is ipi(R), which satisfies 
TLipi = Eiipi, is expanded in bulk Bloch states u„ ^ O^-) °f 
material M. [It should be noted that the bulk materials 
(M) can be strainediSi] Namely, 



To calculate the energies and wavefunctions of electron 
and hole states in the quatum dot, we use the empir- 
ical pseudotential method of Wang and Zungeri2i This 
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Hydrostatic distortion and biaxial strain 
Lens-shaped Ir^ 6 Ga 04 As/GaAs; b=252A 



h=35A 



h=50A 




FIG. 2: Hydrostatic distortion [Tr(e); top] and biaxial strain (bottom) as a function of dot's height plotted on a planar section 
that is normal to [100] and bisects the dot (see Fig. 0. Positions are measured in units of aGaAs = 5.65 A. The wetting layer 
(WL) strain appears prominently. Alloy fluctuations make the strain profiles sample dependent. Two robust features emerge: 
(i) The biaxial strain at the top of the dot increases with height and slowly decays off. The taller the dot the faster this decay. 
Tr(e) shows a similar trend, (ii) The underlying zincblende atomic structure causes the strain to rapidly drop to zero along the 
[Oil] and [Oil] directions. Arrows indicate this feature. 



where, n and k indicate the band-index and wave-vector 
of the Bloch state, respectively; and N is the number of 
primary cells contained in a supercell that encloses the 
quantum dot. Thus, by diagonalizing the Hamiltonian 
[Eq. ©] in the Bloch states basis we find the coefficients 

^M-nk- The calculated wavefunctions ^j(R) are 2-fold, 
Kramers degenerate, so we have omitted the spin index 
a. This representation for the wave function is very dif- 
ferent from the familiar k-p method in that in the latter 
approach the basis set is constructed only from states 
near the Brillouin zone center (r), while here we use a 
full-zone description. Further, in k-p one is restricted to 
just one or two Bloch bands at T (8 x 8 representating 
two Bloch bands) while here we consider n Bloch bands. 
In Appendix^ we present an assessment of the conver- 
gence of energy levels £i as a function of the expansion 
parameters in Eq. JUJl. 

D. Strain-modified band offsets 

While we describe the effects of strain on the elec- 
tronic structure atomisticallyas shown in Eq. |3J here, 
for illustrative purposes only, we calculate the conduc- 
tion (electron) and valence (hole) strain-modified band 
offsets (confining potentials) present in the quantum dot 
by coupling strain to kp-like equations. At each 8-atom 
unit cell in the supercell we diagonalize the conduction 
and valence (including spin-orbit coupling) band strain 



Hamiltonians that Wei and Zunger put forward in Rcf. 
IlH Figure shows the calculated strain-modified band 

offsets for electrons and holes along the line OO' indi- 
cated in Fig. 2] Electron band offset appears in the up- 
per panel, whereas thick and thin lines in the lower panel 
show the first and second holes offsets, respectively. We 
note the following features, (i) The band offsets (Fig. 
|3J) inherit the jagged nature of the strain fields (Fig. 0). 
Inside the quantum dot (region "D" in Fig. , alloy fluc- 
tuations lead to a small mixing in the heavy-hole (HH) 
and light-hole (LH) character of the band offsets, (ii) Re- 
gardless of the dot height, the higher energy hole band 
offset has HH character inside the dot and LH outside. 
Conversely, the lower energy hole band offset has LH 
character inside and HH outside, (hi) The increase of 
biaxial strain at the dot-matrix interface that occurs as 
the height increases (Fig. [2| is reflected in the increase 
(decrease) of the higher (lower) energy hole band offset. 
(See arrows in Fig. |3J|. In particular, for tall dots, the 
decrease of the lower energy hole offset at the interface 
leads to the formation of a pocket in the band offset. 



E. Orbital and Bloch character of wavefunctions 

The single-band effective mass model predicts that in 
the continuum limit, in which lens-shaped dots have full 
axial symmetry, the confined energy levels will form shells 
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Strain-modified band offsets 
Lens-shaped ln 6 Ga 04 As/GaAs; b=252A 
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FIG. 3: Strain-modified band offsets for electrons (upper panel, e), and first and second holes (lower panel; thick and thin 

lines, respectively) along the OO' || [001] as a function of height. (See Fig. In the dot region (indicated with a D) the 
offsets are irregular (jagged) due to alloy fluctuations. These irregularities are more noticeable as the height increases. Outside 
D the heavy-hole (HH) and light-hole (LH) character reverses. Arrows indicate "pocket" formation in the second hole band 
offset as the dot's height increases. 



of degenerate states with definite angular momentum 
I = 0, 1, 2, • • -. 44 Thus, the orbital character of the con- 
fined levels will be pure "s" , "p" , "d" , etc. However, when 
the correct symmetry of the dot is taken in to account, 
as it is the case when using our atomistic pseudopoten- 
tial method, the confined energy levels can have mixed 
orbital character. To quantify the degree of mixing, we 
analyse the single-particle wavefunctions t/>j(R) by pro- 
jecting their envelopes on eigenfunctions of the axial (\\z, 
see Fig. angular momentum e jm< ^/\/27r and determine 
the orbital (angular) character. The latter given by 



The Bloch character of the hole wavefunctions t^i (R) 
in the quantum dot depends on strain as well as band cou- 
pling. To study the extent to which the Bloch character 
is HH, LH, and SO, we proceed similarly and project the 
envelope of tp^ (R) on the total-angular-momentum basis 
| J, J z ) ({|3/2, ±3/2), |3/2, ±1/2), 1 1/2, ±1/2)}). 



F. Exciton energy levels are obtained via the 
screened configuration interaction approach, not 
perturbation theory 



dz I dpp 
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We calculate the exciton energy levels E v by using the 
configuration interaction method as proposed in Ref. l42l 
Briefly, in this method, the exciton states 



In Eq. Q, we have used cylindrical coordinates, i. e. 
R = {p,4>,z), and written the T-derived envelope func- 
tion as 
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where {\eihj)} denotes a basis of single-substitution 
Slater determinants (configurations), in which an elec- 



tron is promoted from if)j n> (R.) to if)^ e '(R.). The CI 
method would deliver the exact exciton ground and ex- 
cited states in the case of a complete basis. However, we 
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truncate this basis and consider all the possible config- 
urations build out of n e electron and hole confined 
states. The coefficients Cy arise from the diagonaliza- 
tion of the exciton Hamiltonian. The direct (Jy = T^) 



and exchange (Kij = r?'^) electron-hole Coulomb in- 
tegrals that enter the calculation are derived from the 
Coulomb scattering matrix elements 



1 k,i — 



dRdR' 



e(R, R')|R- R'| 



(10) 



r 



The microscopic, phenomenological dielectric constant 
e(R, R') that screens the interaction is calculated within 
the Thomas-Fermi model proposed by Resta4£ We do 
not use simple perturbation theory where the exciton is 
described via a Coulomb-corrected single-particle band 
gap (£ I — £j) — Jij: here, £f and £^ are the energies of 
electron level i and hole level j, respectively. 



III. RESULTS 

We now present the effects of height and composition 
on the energies of confined levels and their splittings, 
wavefunctions of selected confined electron and hole lev- 
els, and the lowest transition energy of the exciton. 



A. Energies of confined levels: "p" levels split even 
for ideal lens-shape dots 

The confined electron and hole energy levels are re- 
spectively labeled as £f and where i — 0, 1, • • • is 
an orbital quantum number. The corresponding wave- 
functions are il)n\r) and ip^fa). Each of the con- 
fined states are 2-fold, Kramers degenerate. We la- 
bel the i = electron and hole states as LUMO and 
HOMO, respectively. The first 20 electron and hole en- 
ergy levels of an Ino.6Gao.4As/GaAs dot appear in Fig- 
ure 0] The electron and hole energies are measured, 
respectively, from the bottom of the conduction band 
(CBM) E c (GaAs) = -4.093 eV (calculated bulk electron 
affinity) and from the top of the valence band (VBM) 
E v (GaAs) — — 5.620 eV (calculated bulk ionization po- 
tential) of bulk GaAs. For all heights, the confined 
electron states form groups (shells) of non-degenerate 
levels. 44 In turn, this shell structure can only be identi- 
fied for the first few hole levels (near HOMO) in shallow 
dots (up to h — 50 A). For taller dots, the holes show no 
shell structure. The number of confined electron states 
is significantly smaller than that of holes. While for the 
tallest dot there are 10 confined electron states, more 
than 150 hole states are confined in all the considered 
dots. When comparing the results for h = 35 A with the 
first 20 electron and hole energy levels (not shown) in a 



lens-shaped, non-alloyed InAs/GaAs dot, we find that (i) 
the electron levels form shells that have a bigger average 
separation (55meV vs 45meV in Ino.6Gao.4As/GaAs), 
(ii) the hole energy level structure near the HOMO is 
significantly different (see splittings in Fig. |5J) and (iii) 
the single-particle gap £q — £g is smaller (see Fig. ITUfr . 

Energy splittings. — The "s-p" energy splitting £\ — £q 
and "p-p" splitting £ 2 — £\ for electrons and holes are 
shown in Figure [7| as a function of the dot height. (The 
dots' single-particle gap is also indicated.) Two features 
emerge. First, "s-p" splitting for electrons is bigger than 
for holes. For electrons, the magnitude of the splitting 
is about 45meV and depends weakly on height. On 
the other hand, for holes, this splitting changes from 
~ 20 meV to nearly 5 meV when height changes from 
20 A to 75 A. Second, the "p-p" splitting shows the op- 
posite behavior. Namely, for holes the magnitude of 
the splitting remains nearly constant at approximately 
4meV, and for electrons it changes non-monotonically. 
It should be noted that the electronic "p-p" splitting is 
sensitive to alloy fluctuations and it can change by al- 
most a factor of two by changing the alloy realization in 
the dotM. 

For comparison, "s-p" and "p-p" splittings in a lens- 
shaped, non-alloyed InAs/GaAs dots with b = 252 A also 
appear in Fig. \5\ The pure dot has a "s-p" splitting 
for holes that is nearly the same as in the h = 35 A 
Ino.6Gao.4As/GaAs dot, while this splitting for electrons 
is nearly 20% bigger and depends weakly on height. In 
contrast, the "p-p" splitting for holes in the pure (non- 
alloyed InAs/GaAs) dot is about twice as big as in the 
h = 35 A alloy dot, and for electrons the "p-p" splitting 
depends weakly on height and is similar in magnitude 
to the splittings in alloy dots. It should be noted that 
in pure (non-alloyed InAs/GaAs) dots the hole energy 
levels undergo a localization crossover (see discussion in 
Sec. IIII CJ|) for tall dots, which render meaningless the 
notion of "s-p" and "p-p" splittings. (Such a localization 
crossover is absent for electrons.) For this reason, we have 
compared the splittings for the holes in a shallow (h = 
35 A) non-alloyed InAs / GaAs dot only, while comparing 
the splitting for electrons in a range of heights. 
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Electron and hole energy levels 
Lens-shaped ln 6 Ga 4 As/GaAs; b=252A 
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FIG. 4: First 20 electron (top panel) and hole (bottom) en- 
ergy levels Sn and £„ for Ino.6Gao.4As/GaAs dots with dif- 
ferent heights. The electron and hole energies are measured 
from GaAs CBM [E c {GaAs) = -4.093 eV] and GaAs VBM 
[E v (GaAs) = — 5.620 eV], respectively. As height increases, 
the number of confined levels (£„ < and > 0) increases 
and the single-particle gap (£q —£q) decreases. The confined 
electronic energy levels group in non-degenerate shells for all 
dots, "s", "p", and "d" indicate the predominant orbital char- 
acter of selected states. For holes, up to h = 50 A the second 
and third levels are "p"-like, while for larger heights these 
levels are "s"-"p" hybridized. 



B. Wavefunctions of confined states in 
Ino.6Gao.4 As/GaAs and non-alloyed InAs/GaAs 
dots: mixed "s" + "p" character and "HH" + "LH" 
character 



Figure[5]compares the wavefunctions of LUMO (lowest 
electron) and HOMO (highest hole) states as a function 
of height in Inn.eGan. 4 As/GaAs dots with base b — 252 A. 
To make the comparison, we plot isosurfaces that en- 
close 75% of the total charge density, and show contour 
plots taken at lnm above the dot's base. In addition, 
the V -orbital character of the LUMO and HOMO is 
indicated as well as the HOMO's heavy- and light-hole 
character. Further, for each height, the energy £q of 
LUMO is shown relative to E c (GaAs) and the energy of 
HOMO £q relative to E v (GaAs). Prominent results are 
the following. LUMO: The lateral spatial extent of the 



"s-p" and "p-p" electron and hole energy splittings 
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FIG. 5: Height dependence of the "s-p" [£\ — £0) and "p- 
p" (£2 — £1) energy splittings for electron and hole levels in 
Ino.6Gao.4As/GaAs. The single-particle gap is also shown. 
The "s-p" splitting for electrons is bigger than for holes; in 
turn, the "p-p" splitting shows the opposite trend. The height 
dependence of the "p-p" splitting for electrons is not mono- 
tonic, due to alloy fluctuations, while for holes this splitting 
remains nearly unchanged. For comparison, we show the "s- 
p" and "p-p" splittings for electron (solid squares) and hole 
(open squares) in non-alloyed InAs/GaAs dots. 



wavefunction depends weakly on height. In contrast, as 
height decreases, the wavefunction extends into the bar- 
rier along the vertical direction. The "s" -orbital char- 
acter of LUMO remains at nearly 90%. HOMO: Wave- 
functions are more sensitive to height, showing a spatial 
extension that gets reduced significantly both in the lat- 
eral and vertical direction when height changes from 75 A 
to 20 A. This reduction leads to a strong localization at 
the center of the dot for h — 20 A. The "s" -orbital charac- 
ter of the HOMO remains at about 90% up to h = 50 A, 
for taller dots the character of the HOMO becomes "s" 
and "p" mixed. This mixing reflects the reduction of 
the hole charge density near the center of the dot. The 
heavy- (HH) and light-hole (LH) character also change 
with height in a similar manner as the "s" and "p" char- 
acter. Namely, for the three smaller dots (20 A, 35 A, 
50 A) the HOMO is mostly of the HH type, but as the 
height increases the LH character increases. We expect 
this behavior since the LH band-offset increases within 
the dot as height increases. (See Fig. 0) 

Excited hole states. — Figure [7| compares wavefunctions 
for excited hole states; namely, second (HOMO— 1) and 
third (HOMO-2) hole states. As in Fig. we show 
isosurfaces and contour plots, present the orbital and 
HH/LH character of these states, and the energies £ f and 
£\ for HOMO-1 and HOMO-2, respectively. The char- 
acter of HOMO— 1 and HOMO— 2 are nearly the same 
at all heights. In addition, these states have a domi- 
nant "p" -orbital character for all heights, regardless the 
absence of "p" -shell structure for tall dots (h = 65 A 
and 75 A). (See Fig. As in the case of the HOMO 
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FIG. 6: Height dependence of lowest electron level (LUMO) 
and highest hole level (HOMO) wavefunctions. Isosurfaces en- 
close 75% of the total charge density. Countours are taken at a 
plane lnm above the base. The energy relative to E c (GaAs) 
and E v (GaAs) is shown. The "s-" and "p-"orbital charac- 
ter (see Eq. |7J , and the heavy-hole (HH) and light- hole (HH) 
character are indicated. For the tallest dot, the wavefunctions 
are entirely confined within the dot and the effective size of 
the dot — spatial extent of the wavefunctions — is nearly the 
same for LUMO and HOMO. For shorter dots the effective 
size gets significantly reduced for the HOMO while it remains 
nearly unchanged for the LUMO. 



state, the HOMO— 1 and HOMO— 2 states have increas- 
ing light-hole character with increasing height. However, 
the percentage of light-hole character is almost twice that 
of HOMO. For instance, at h = 75 A, the LH character 
of the HOMO-2 is 17% while the HOMO LH character 
is 11%. 
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FIG. 7: Wavefunctions of the first 2 excited hole levels 
HOMO-1 and HOMO-2 for different dot's height. As in Fig. 
|S| the isosurface encloses 75% of the total charge density and 
the countours are taken at lnm above the base. Labels in- 
dicate the "p" orbital character, the heavy- and light-hole 
characters and the energy of the states relative to E v (GaAs). 



C. Strain-driven hole localization 

For dots containing both Ga and In 
(Ino.6Gao.4As/GaAs) with b = 252 A and heights 
in the rang e 20-75A, we have shown that the wavefunc- 
tions of HOMO and HOMO-1 as well as 
other low-lying excited hole states, are localized inside 
the dot (see Figs. El and 0) and that these levels have 
a sizeable "s-p" splitting E\ — £q (see Fig. |5j. In 
contrast, in lens-shaped pure, non-alloyed InAs/GaAs 
quantum dots, localization of the low-lying hole states 
at the dot interface develops as the height of the dot 
increases. In addition, HOMO and HOMO— 1 become 
nearly degenerate [e. g., — E v (GaAs) = 256.2 meV 
and E\ - E v (GaAs) = 255.7meV at h = 75A], as 
well as HOMO-2 and HOMO-3. (It should also be 
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HOMO and HOMO-1 wavefunctions 
Lens-shaped InAs/GaAs; b=252A 




FIG. 8: HOMO and HOMO-1 wavefunctions in pure, non- 
alloyed InAs/GaAs lens-shaped quantum dot as a function of 
height. The dots have the same base diamater b — 252A. 
The aspect ratio b/h is shown, as well as the energy relative 
to E v (GaAs). As before, the isosurfaces enclose 75% of the 
charge density and the countours are taken at 1 nm above the 
base. As height increases, interface localization takes plave, 
and the HOMO and HOMO— 1 become nearly degenerate. 



offset) along OO' and that the band offsets values de- 
crease. Further, the band offset become assymetric with 
the values at the top of the dot smaller that at the base. 
As a consequence, it becomes clear why the energy of 
the HOMO and HOMO-1 become bigger and we also 
expect the wavefunctions of these states to be localized 
near the base of the dot. It is important to mention 
that in the calculation of the energy of the HOMO and 
HOMO— 1 we do not utilize the band offset we present 
in Fig. EJa). Second, the band offsets shown in Fig. |5| 
present a "crown" structure at the dot-matrix interfaced 
This crown structure becomes more significant as height 
increases, as a consequence of the values of band offset 
along the [110] and [110] direction presenting a weak de- 
pendence on height. Finally, the localization of the low- 
lying hole wave functions that develops as the height of 
pure, non-alloyed InAs/GaAs lens-shaped quantum dots 
is a results of the peculiar characteristics of the higher en- 
ergy valence-band offsets, which in turn are determined 
by the biaxial strain profile in the nanostructure. 

Finally, it should be noticed that electron states do 
not experience interface localization; instead, they con- 
tinue to form shells and the levels in these shells have a 
predominant, not-mixed orbital character ("s", "p", • • •). 
The reason for the absence of electron localization is that 
electron states are sensitive to the (hydrostatic) isotropic 
strain rather than to the biaxial strain. 



noted that for the flat dot the energies of HOMO and 
HOMO— 1 are bigger than in the alloy dot of the same 
size.) Remarkably, HOMO and HOMO-1 are polarized 
along [110] (Fig. |SJ while HOMO-2 and HOMO-3 
(not shown) are polarized along [110]. Figure |H1 shows 
the development of this interfacial localization and the 
energy of HOMO and HOMO-1 relative to E v (GaAs). 
As in Figs. and isosurfaces enclose 75% of the 
HOMO and HOMO-1 total charge density. 

It should be noted that the pure, non-alloyed 
InAs/GaAs dots have the same base size and same height 
range as the alloy dots we have previously discussed. 
Hence, hole localization does not have its origin on geo- 
metrical aspects of the quantum dots. Instead, the local- 
ization of the hole wavefunctions is driven by the biaxial 
strain present in the nanostructure (QD+GaAs matrix). 
To explain this result, we plot in Fig. E{a) the valence 
band offsets of heavy- hole (HH) and light-hole (LH) char- 
acter (inside the dot) along the line OO' (see Figs. Q] 
and El) and the energy of HOMO and HOMO-1 (thick 
dashes). In addition, we plot in Fig. Etb) a 3-dimensional 
rendering of the higher energy hole band offset values at 

a plane normal to OO' , cut slightly above the base of 
the dot. First, we see that as the dot heigh increases 
the pocket structure [indicated with arrows in Fig. EUa)] 
that appears in the higher energy band offset dramati- 
cally widens the confining potential (given by the band 



D. Exciton transitions 

Figure ITU1 shows the lowest bright transition (gap) of 
the exciton as a function of dot height, calculated at the 
single-particle (SP) level (£q—£q) and by using the many- 
body, configuration interaction (CI) method. 42 In the lat- 
ter, we use a basis consisting of 12 electron and 20 hole 
confined levels (n e — 12, nh — 20) for heights h = 35- 
75 A and (n e = 6, rih = 20) for h = 20 A. As expected, 
the transition energy decreases as height increases, due to 
confinement. However, the scaling with height differs sig- 
nificantly from predictions of single-band, effective mass 
(~ /i~ 2 )4£ Namely, by fitting our results for the height 
dependence of the gap to the function a + b/h 1 we find 
■y SP = 0.95 and j CI = 1.09. The value of the gap for 
large heights correspond, respectively, to a SP — 1.117eV 
and a CI — 1.116 eV in the single-particle and CI ap- 
proach. We expect the SP and CI scaling exponents to 
be different, as in the single-particle calculation the scal- 
ing is dictated by the scaling of the LUMO and HOMO, 
whereas in the CI calculation the electron-hole matrix 
elements are also included and the magnitude of these 
matrix elements decreases with height. The discrepancy 
between the scaling (7) in single-band effective mass and 
both our SP and CI calculations can be attributed to non- 
trivial effects that are naturally accounted for within our 
atomistic approach such as non-parabolicity and multi- 
band effects, and the position-dependent strain present in 
the nanostructure. The values of a at large heights is also 
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FIG. 9: Strain-modified band offsets in a pure, lens-shaped, non-alloyed InAs/GaAs quantum dot along the OO' direction (a) 
and in a plane normal to latter cut at a height slightly above the base of the dot (b). Dashes in panel (a) correspond to the 
HOMO and HOMO — 1 energy levels. The splitting among these levels decreases with the dot height. Arrows indicate pocket 
formation (see also Fig. in the band offset due to strain accumulation at the dot-matrix interface at the top of the dot. 



expected to be different, and the difference is attribute 
to correlation effects that are present in the many-bod^ 
configuration approach. Thus, we expect a small diffci 
ence (~ meV) between a and a CI . 

For comparison, we also present results for the lowes 
exciton transition [squares, open (SP) and solid (CI)] in 
lens-shaped, non-alloyed InAs/GaAs quantum dot. Th 
values of the transition energies are smaller than in th 
alloy dots. 



IV. SUMMARY 

By using a high-level atomistic approach, we hav 
predicted spectroscopic characteristics of self- assemble 
Ini_ x Ga x As/GaAs quantum dots as a function of heigl 
and composition. Several prominent features emerged. 

(i) The biaxial strain at the dot-GaAs matrix intei 
face increases with height, whereas, regardless of heighi 
the strain is negligibly small along the [011], [011] an 
crystallographically equivalent directions. 

(ii) Regardless of height and composition, the confrncu 
electron energy levels group in shells of nearly degen- 
erate states. The average energy splitting among these 
shells depends weakly on height; however, this splitting is 
larger in pure, non-alloyed InAs/GaAs. In contrast, the 
confined hole energy levels form shells only in flat dots 
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FIG. 10: Exciton gap as a function of height. Single-particle 
(open circles) and configuration interaction (CI, solid circled) 
results are shown. The CI basis is (n e = 12, n/, = 20) for 
h = 35-75 A and (n e = 6, n h = 20) for h = 20 A. Dashed lines 
represent fits to the function a + b/h 1 . Fitting parameters are 
indicated. For comparison, we present results [squares; open 
(SP) and solid (CI)] for a pure, non-alloyed InAs/GaAs dot. 
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and near the highest hole level (HOMO). 

(iii) In alloy dots, the electrons' 1 "s-p" splitting de- 
pends weakly on height, while the "p-p" splitting de- 
pends non-monotonically — due to alloy fluctuations. In 
pure, non-alloyed InAs/GaAs dots, both these splittings 
depend weakly on height. Further, the "s-p" splitting 
is larger while the "p-p" has nearly the same magnitude. 
For holes levels in alloy dots, the "s-p" splitting decreases 
with increasing height (the splitting in tall dots being 
about 4 times smaller than in flat dots), whereas the 
"p-p" splitting remains nearly unchaged. Shallow pure 
dots have a "s-p" splitting of nearly the same magnitude, 
whereas the "p-p" splitting is about three times larger. 

(iv) As height increases, the "s" and "p" character of 
the wavefunction of the HOMO becomes mixed, and so 
does the heavy- and light-hole character. 

(v) In alloy dots, regardless of height, the wavefunc- 
tion of low- lying (near the HOMO) hole states are lo- 
calized inside the dot. Remarkably, in pure, non-alloyed 
InAs/GaAs dots, as the dot height increases, these states 
become localized at the dot-matrix interface and nearly 
degenerate. Further, the localized states are polarized 
along [110] (HOMO and HOMO-1) and [110] (HOMO-2 
and HOMO-3). This localization effect is driven by the 
peculiarities of the biaxial strain present in the nanos- 
tructure. 

(vi) The lowest exciton transition energy (gap) de- 
creases with height, but the scaling (roughly ~ /i _1 ) dif- 
fers significantly from the prediction of single-band effec- 
tive mass (~ h~ 2 ). 

The study we presented here may be used to bridge 
spectroscopy results with theory without the need for 
severe theoretical approximations. 
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APPENDIX A: CONVERGENCE 

The LCBB method contains several convergence pa- 
rameters: supercell size, number of bulk materials (M), 
number of bands (n), and number of k points. The choice 
of bulk materials and number of bands can be physically 
motivated and justified, whereas the choice of supercell 
size and number of k points is not as clear. For this 
reason, we conducted a converge assesment on supercell 
size and number of k points. As for the bulk materials, 
we choose (i) unstrained GaAs and (ii) InAs subject to 
strain values of e xx = e yy = —0.06, and e zz = +0.04. 



These values are typical of strain inside a non-alloyed 
InAs/GaAs lens-shaped quantum dot. For the bands n, 
when simulating the electron and hole energy levels we 
use respectively the lowest conduction band and the three 
highest valence bands of both bulk materials. 

Supercell size. — We consider cubic supercells of size 
L x xL y xL Zl where L t (t = x,y,z) is measured in units 
of the lattice parameter of bulk GaAs aaaAs- Based on 
strain calculations (see above), we selected a supercell 
size such that the strain profile within the quantum dot 
remained unchanged upon changing size. Notwithstand- 
ing, it should be noted that although the strain may be 
converged within the dot the strain values at the bound- 
ary of the supercell may not be zero due to the slow 
relaxation of the strain fields. This residual strain intro- 
duces small changes (~ meV) in the energy levels. For 
instance, Figure ^shows the effect of reducing the resid- 
ual strain by changing the supercell size from 72x72x40 
to 72x72x48 for a h = 20 A dot. 

Number of k points. — The k-point meshes enclose all 
the k vectors around the T point that lie within an ellip- 
soid with main axis equal to (2ir/L x ) P, (2n/ L y ) Q, and 
(2tt/L z ) K, respectively. P, Q, and K give the number of 
k points taken along each cartesian direction. In our con- 
vergence assesment P = 6, 8, Q = P, and K ranges from 
8-18. Figure ITT1 shows the electron energy levels for sev- 
eral dots as a function of k-point mesh. We distinguish 
several features, (a) When compared with the energy of 
the levels, differences in energies are quicker to converge, 
(b) High-energy levels require larger k-point meshes to 
converge within a given threshold, (c) The taller the 
quantum dots the smaller the number of k points needed 
to converge, (d) The use of (6, 6, K) and (8, 8, K) k- 
point meshes results in similar energies for the lowest 
electron level. The discrepancy between the predictions 
made with these two meshes increases for high-energy 
levels. 

To present details of the convergence with respect to 
k-point mesh, Figure El shows a convergence assessment 
of HOMO and LUMO's energy in flat (h = 20 A) and tall 
(h = 75 A) dots. The energy of HOMO converges visibly 
quicker than that of LUMO. This behavior holds for all 
hole energy levels, and arises from the higher size (three 
Fi5„ bands) of the basis used in the simulation. 

To summarize, the results we present in this work for 
each quantum dot derive from simulations with k-point 
meshes that provide with energy levels converged within 
2meV. Table [I] shows the supercell sizes and k-point 
meshes we have used to simulate the electronic struc- 
ture of each of the dots we have studied. It should be 
noted that energy differences are converged to much less 
than this lower bound. 
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TABLE I: Summary of supercell sizes L x xL y xL z (lengths in 
units of aaaAs) and k-point meshes (Q, Q, K) used in the 
simulations. Dot's base b = 252 A. (c.f. Fig. for dot's 
geometry.) 
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